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Abstract 

The Beth Uhlenbeck formula gives an exact (quantum) expression of 
the second virial correction to the equation of state of a dilute gas. We show 
how this result can be extended to arbitrary degeneracy. For this purpose 
we develop a formalism based on the use of Ursell operators, which contain 
no symmetrization in themselves (they can be used for both a system of 
identical particles or an auxiliary system of distinguishable particles). A 
concise expression generalizing the Beth Uhlenbeck formula is obtained, 
which is equally valid for bosons and fermions. The formalism is not es- 
sentially limited to dilute gases, but can also be applied to the study of 
denser systems, for instance those undergoing a Bose Einstein transition; 
this will be the subject of a forthcoming article. 



1 Introduction 



The Beth Uhlenbeck formula [JTJ] gives an exact expression of the second virial 
coefficient as a function of all collision phase shifts associated with the interaction 
potential. This is a remarkable result which, in a concise formula including no 
phenomenological constant, relates macroscopic and microscopic quantities: on 
one hand the density, the temperature, and the pressure of a gas, on the other 
the quantum phase shifts, which are exactly calculable from the Schrodinger 
equation. A closer inspection shows that the pressure correction is the sum of two 
contributions, one of pure statistics, and a second arising from the interactions; 
the former is actually nothing but the first order term of the expansion in n\ 3 
of the ideal gas equation of state. It is therefore not surprising that the validity 
of the formula should require two conditions, that the gas be dilute with respect 
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to both degeneracy effects (n\ 3 <C 1) and interaction effects as well (nb 3 <C 1, 
where b is the range of the interaction potential). A natural question is whether 
one can release the first of these conditions and study gaseous systems over a 
larger range, going continuously from the classical region where nA 3 <C 1 to the 
quantum region where this quantity is comparable to one. One would just assume 
that the gas remains always dilute in terms of the interactions. As a matter of 
fact, most classical textbooks on quantum statistical mechanics include a study 
of gaseous systems that are degenerate and imperfect. The methods they use 
are based on various forms of perturbation theory; for instance they introduce 
an auxiliary pseudopotential which replaces the real interaction potential 0, or 
they make a direct substitution of a collision T matrix for this potential @]. 
This well known approaches are, nevertheless, different from that of Beth and 
Uhlenbeck whose calculation is exact (of course they limit the calculation to the 
first density correction to the equation of state, but this is precisely the definition 
of the second virial correction). In essence, the cluster expansion method they 
use is a method providing directly density expansions, as opposed to interaction 
expansions, which require at least partial resummations of diagrams in powers of 
the interaction potential || to reconstruct the T matrix from V and to obtain 
density expansions. Actually, the exact treatment of the relative motion of two 
particles is explicit in the Beth Uhlenbeck formula, which contains a thermal 
exponential of the two particle hamiltonian ensuring that thermal equilibrium has 
been reached, in other words that the binary correlations between the particles 
are treated properly. Technically, a characteristic of this approach is that no 
expression which becomes infinite for repulsive hard cores is ever written, even 
at some intermediate stage. It therefore seems to be an attractive possibility to 
try and extend the method to degenerate systems, if only to compare the result 
of both methods, interaction and density expansions. This is the subject of the 
present article; we show that concise and exact formulas can indeed be written 
that provide a natural generalization of the Beth Uhlenbeck formula, without 
requiring complicated algebra. 

For this purpose, we will make use of a generalization of the quantum UrsellQ 
cluster functions introduced by Kahn and Uhlenbeck 0: Ursell operators which 
intrinsically do not contain statistics (as opposed to the usual quantum Ursell 
functions H[f§, which are already symmetrized). The action of the operators is 
defined, not only in the space of symmetric or antisymmetric states of the real 
system, but also in the larger space of an auxiliary system of distinguishable par- 
ticles. If the system is dilute in terms of the interactions but not of statistics, 
this makes it possible to limit the calculation to low order Ursell operators while, 
with the usual (fully symmetrized) Ursell functions, one would need to include 

1 The initial introduction of the Ursell functions was made by him within classical statistical 
mechanics R], and generalized ten years later to quantum mechanics by Kahn and Uhlenbeck. 
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higher and higher orders with increasing degeneracy^. Of course, this method 
implies that we have to give up well-known advantages of the formalism of sec- 
ond quantization, but this is the price to pay for the treatment of interactions 
and statistics in completely independent steps. It also means that an explicit 
symmetrization of the wave functions becomes indispensable at some point and, 
moreover, that no approximation whatsoever can be made at this step: otherwise 
the possibility of treating strongly degenerate systems would be lost. Fortunately 
it turns out that the symmetrization operation can indeed be performed exactly. 
This is done by introducing simple products of operators, corresponding to ex- 
change cycles, or more generally simple functions of operators (fractions) that 
correspond to summations over the size of these cycles up to infinity. Physically, 
this emphasizes the relation which exists between the properties of a system, its 
pressure for instance, and the size of exchange cycles [10] that take place in it. 

Another feature of the method is that it naturally provides expressions for 
the one and two body density operators, which is of course convenient if one is 
interested in a detailed study of correlations. The formalism is useful beyond the 
evaluation of virial corrections only; this is why we develop it in what follows 
beyond what is strictly necessary to generalize the Beth Uhlenbeck formula. A 
preliminary report on this work has already been given in flOf ; applications to 



denser systems will be the object of a forthcoming article |11 



2 Ursell operators 

In the canonical ensemble, the partition function of a system of indistin- 
guishable particles is given by a trace inside the space of symmetrical (or anti- 
symmetrical) states. Nevertheless, with the help of the projectors S and A onto 
these (sub) spaces, the trace can be extended to a larger space, which is associated 
with the same number of distinguishable particles: 

Z N = Tr{K N S A } (1) 

In this equation, S applies for bosons while A applies for fermions, and the 
operators K N are defined by: 



K x = exp-/?# (l) 

K 2 = exp-(3[H (l) + H (2)+V 12 ] 

K N = exp -p [E£i#o(i) + £i>,%- 



(2) 



2 As discussed in p0[ , when a Bose Einstein condensation takes place, the Ursell operators 
have no singular variation, which is another way to see that they are not sensitive to degeneracy 
by themselves. 



3 



with usual notation: Hq(1) is the one particle hamiltonian including kinetic en- 
ergy and an external potential (if there is any in the problem), is the in- 
teraction potential between particles with labels i and j. We now use cluster 
techniques to expand operators exactly in the same way as one usually does for 
functions. The Ursell operators U\ are therefore defined according to: 



Ui = K x 

C/ 2 (l,2) = K 2 (l, 2) -7^(1)^(2) 
E7 3 (l,2,3) = K 3 {1, 2, 3) -#2(1,2)7^(3) - K 2 (2, 3)^(1) 
-K 2 (3, 1)^(2) + 2K 1 (l)K 2 (2)K l (3) 



(3) 



etc. (I = 1,2, ..N). Conversely, in terms of the Ursell operators, the N particle 
operator Km can be written in the form: 

Kn= E E Ui(-)Ui( - )-Ui(-) x U*(; W- , -)-^(-, ■) x£/ 3 (., ., ■) ■■■■ (4) 

where the first summation is made on all possible ways to decompose the number 
of particles as: 

N = E lm 'i (5) 
i 

The second summation correspond to all non-equivalent^ ways to distribute the 
iV particles into the variables of the Ursell operators, symbolized by dots in @. 
It is convenient to simplify the two summations into one: 

E E ^E (6) 

{U} 

which is meant over all non-equivalent ways to distribute the particles into various 
sequences of U's. There is little difference between our definitions and those of 
section 4.2 of [|J or ||: we use operators instead of symmetrized functions and, 
more importantly, the action of these operators are defined, not only within the 
state space that is appropriate for bosons or fermions, but also in the larger space 
obtained by tensor product that occurs for distinguishable particles. Hence the 
need for an explicit inclusion of S or A in (fl|). 

We now decompose these operators into a sum of permutations P a that, in 



3 For instance 1/2(1, 2)1/3(3, 4, 5) and 1/2(2, 1)173(4, 3, 5) correspond to equivalent distributions 
of 5 particles, since the very definition of the Ursell operators implies that the order inside each 
Uk is irrelevant. 
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turn, we decompose into independent cycles C of particles^} 

! = 4 E E ^ C 2 (., .)C 2 (., .)C a (, .) xC 8 (, ., .)•••• (7) 

A ^'| m |{D) s « ' s ' 

mi factors rri2 factors 

where the first summation is similar to that of @, while the second corresponds 
to all non-equivalentQ ways to distribute numbers ranging from one to N into the 
variables of the various C"s. In this equation, the operator S applies for bosons 
with i] = +1, while A with 77 = — 1 applies for fermions. We also simplify the 
notation into: 

E E =► E (8) 

{mi}{D} {P a } 

We can now insert (f|) and (0) into (0) and obtain, within a double summation, 
numbers that are traces calculated in the space of distinguishable particles, i.e. 
in the ordinary tensor product of iV single particle state spaces. Inside most 
of the terms of the summation, factorization into traces taken inside smaller 
subspaces occur. For instance, if the term in question contains particle number % 
contained at the same time in a U\ operator as well as in a C%, the contribution 
of that particle completely separates by introducing the simple number Tr{U\}. 
Or, if n particles are all in the same U n but all in separate C\S, this group 
of particles contributes by a factor Tri„ n {U n } ; if they are all in different UxS 
but also contained in one single large cycle C n , their contribution also factorizes 
separately. More generally, in each term of the double summation, particles group 
into clusters (U-C clusters), which associate together all particles that are linked 
either[] by cycles C\ (with I > 1) or by Ursell operators Uf (with / > 1). The 
general term is therefore the product of the contributions of all the clusters that 
it contains and one can write: 

Zw= j^EE II r cluster(i,j,k,..) (9) 
{Pa} {U} clusters 

where k, ..) is the index number of particles contained inside the cluster; the 
number of clusters into which each term of the double summation is factorized 
depends, in general, on this particular term. 

4 We use the same notation as in ref . : Ck(i,j,k,..) denotes a cycle where particle i replaces 
particle j, particle j replaces particle k, etc. This should not be confused with the notation 
P(i,j, k, ...) for a general permutation (not necessarily a cycle) where particle i replaces particle 
1, particle j particle 2, particle k particle 3, etc. Any P can be decomposed into a product of 
C's; for instance, P(2, 1,3) = C 2 (l,2)Ci (3). 

5 Inside every cycle, a circular permutation of the variables has no effect and therefore does 
not affect the permutation P a . 

6 There are therefore two explicitly distinct origins to the clustering of particles in this point 
of view: belonging to the same Uf (with I > 1) or to the same cycle Ci (with / > f). This is 
distinct from usual cluster theories where, either only interactions introduce clustering (as in 
classical statistical mechanics), or the two origins are not explicitely distinguished (as in usual 
quantum cluster theory). 
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3 Diagrams 



Clusters differing only by the numbering of particles that they contain give the 
same contribution. It is therefore useful to reason in terms of diagrams (U-C 
diagrams), which emphasize the way particles are connected through exchange 
cycles and Ursell operators, rather than their numbering. For instance, the first 
diagram will correspond to one particle in a f/i and in a C\ (whatever the num- 
bering of the particle is), and contribute the value Tr{U\} as mentioned above; 
another diagram will introduce the value Tri^-ft^}, etc.; more generally, the 
value of any diagram containing n diag . particles can be written as a trace over the 
variables of n diag . particles numbered arbitrarily: 

IW= ^+^+-Tr 1)2 ,.. ndiag .{C/ 1 (l)x 

x..L/i(n)l7 2 (ri + l,7-i + 2)...C 1 (i)...C 1 (j)C 2 (l, s)...} {LU) 

Here is the number of cycles[] of length k\ the factor r) P2+P4+ ~ arises from 
the factors 77's in (|7|) and corresponds to the contribution of the parity of the 
permutations contained in the diagram to the total permutation of the N particles. 

Now, inside each term of the multiple summation, a given diagram T d i ag may 
occur several times; we then note m d i ag , the number of times it is repeated, and 
we get: 

JV - {P a }{U}diag. 

with the obvious relation: 

N = m diag. X n diag . (12) 

diag. 

Of course identical diagrams appear, not only in the same term of the double 
summation, but also in many different terms. Therefore, if J2f 1 symbolizes 

\ rn diag. J 

a summation over all possible ways to decompose N according to (|P2"D, we can 
also write: 

Z n = ^ E c{m diag .}x J] [r d ia 9 .] mdia3 - (13) 

' {m diag .} dia 9- 

where c {rridiag.} is the numbers of terms in the double summation of ([11]) that 
correspond to this particular decomposition of N. 

To evaluate this number, we have to specify more precisely how the U-C dia- 
grams are constructed graphically. In every cluster, we represent the permutation 
cycles Ck by horizontal lines containing k boxes, or segments, which are available 
to numbered particles. When the corresponding particles are inside Ui's, we do 
not add anything to the diagram; when they are contained inside U^s, we join 



7 Onc obviously has: pi + 2p 2 + ... = ridi 
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the corresponding segments by an additional double line, a triple line for £7 3 's, 
etc. For instance, if we consider the pure exchange cluster: 



r duster . = Tr 1)2 ,.. i7 {C/ 1 (l)?7 1 (2)..C/ 1 (7)C r (l,2,....7)} (14) 

the corresponding diagram will be that of figure 1-a. This kind of linear diagram 
is the only possibility for an ideal gas; to generate a Truster appearing in (f|), it 
must receive a numbered particle in the first segment of the line, in the second 
the particle which is replaced by it under the effect of the permutation cycle, 
in the third the particle which is replaced by that in the second segment, etc. 
Equation ( |i~4|) gives what we will call the "explicit value" of the contribution of 
the diagram; using simple transformations (see ref.JIO] or section |5.1| ), one can 
obtain the simpler "reduced value" Trx |[C/i(l)] 7 |. 

If we now start from a cluster which contains one single V 2 : 

T c i us ter = Tn A 9 {U 2 (l, 2) U x (3) U x (4) U x (5) ..Ux (9)C 3 (1, 3, 4)C 6 (2, 5, ....9)} (15) 

we obtain the diagram shown in figures 1-b; in the same way, the cluster: 

r ' dus ter = Tr lt2i 10 {U 2 (l, 6)U 1 (2)U 1 (3)..U 1 (10)C 10 (1, 2, 3, ., 6, ..10)} (16) 

leads to the diagram of fig. 1-c. These two kinds of diagrams turn out to be 
the only ones that are necessary to generalize the Beth Uhlenbeck formula; we 



give their reduced values in section 5.2. Fig. 1-d shows an example of a diagram 



containing one single C/3, and arising from the trace: 

r = Trx, 2 , ,16 {EMI, 2- 3)L/i(4)EM5)..tfi(13)C 3 (l, 4 > 5 )^(2, 6..12)C 5 (3, ..16)} 

(17) 

and figure 1-e gives another similar example. Clearly, the process can be gener- 
alized to associate a diagram to any cluster, however large and complex. 

But it is not sufficient to construct diagrams, we must also choose explicit 
rules^| ensuring that every T c i uster appearing in where the particles are still 
numbered, will correspond to one single, well defined, diagram; for instance, for 
the cluster written in (p~5f ) , we must decide whether the diagram will be that shown 
in figure 1-b or another where the lowest cycle ("base cycle") is that of length 
6 instead of 3. We will avoid this kind of ambiguity by choosing rules which fix, 
for every cluster associated with a given diagram, where exactly each numbered 
particle should fall into the diagram. These rules ensure that no double counting 
of clusters may occur; they also determine how the geometrical characteristics 
of the diagrams (the lengths of the successive cycles) can be varied, which in 
turn determines the way in which the summations over lengths will be made in 
a second step (next section). For our purpose in this article, which is limited 

8 There is some flexibility in this choice, and here we attempt to take the most convenient 
convention, but it is not necessarily the only possiblity. 
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to a generalization of the Beth Uhlenbeck formula and to the diagrams of figure 
l-a,b and c, we do not need to study the most general case (this discussion is 
given in Appendix I); it will be sufficient to remember that the first particle in 
the lowest cycle ("base cycle") must be, among all particles contained in the Ui 
of highest order I, that which has the lowest index number for the configuration 
to be correct. 

The preceding rules also introduce the "counting factors" corresponding to 
the probability of obtaining a correct representation of a cluster by throwing ran- 
domly numbered particles into a diagram; the counting factors are used below to 
obtain the value of c{mdi ag }- Suppose for instance that we distribute k num- 
bered particles into a linear diagram of the kind shown in fig. 1-a in all possible 
ways; it is clear that there is a probability \jk that the first particles will have 
the lowest index number as required, which leads to the following counting factor 
fdiag. for pure exchange cycles: 

fdiag. — r (18) 

(in other words, if one puts numbered particles into a cycle Cfc(., ., ., ..), one can 
obtain the same permutation k different times). Similarly, it is easy to see that the 
counting factors of the diagrams of figures 1-b and c are 1/2, which corresponds 
to the probability of having the particle numbers contained inside the single U2 
in the correct order. The general value of counting factors is given in Appendix 
I; in many practical situations (when case (ii) of this Appendix does not occur) 
we can ignore the g's in formula (|62|) and use the simpler form: 

fdiag. = [PI M X ImY 1 (19) 

where Im is the largest order of the Ursell operator contained in the diagram 
(one for fig. 1-a, two for figs, b and c, three for figs. 1-d and e) while pi M is the 
number of these Ursell operators in this particular diagram (7 for fig. 1-a, one 
for all the other cases in this figure). Two other examples of diagrams for which 
this formula is valid are given in figure 2; as an illustration, they are explicitly 
calculated in Appendix II. 

We are now in position to calculate c{rridiag.}- This number can be obtained 
by distributing the N particles inside the sites of all the diagrams of a series 
defined by the m's, which can be done in N\ different ways, and counting how 
many times the same term of the double summation of (|TT| ) is obtained. Since 
there are rridiag} ways to interchange the order of all U-C clusters arising from 
the same diagram, there is a first redundancy factor equal to Tiding. ( m diagty that 
comes in. Moreover, there is also only a proportion {fdiag.) mdia3 ' of the obtained 
configurations that is acceptable. Altogether, the net result is: 

c {m diag } = N\ J] [fd ia g] md ^ (20) 

diag. m diag.- 
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4 Two summations 



We now take advantage of the fact that (|20|) contains factorials and, in a second 
step, that the value of fdiag. depends on the topology of the diagrams, but not of 
the size of the chains of Uis that it contains (except for pure exchange sycles). 
This allows us to group together series of terms in (|13"1). The first summation 
is done classically by going to the grand canonical ensemble and defining the 
corresponding partition function by: 

Z 9 . c =Y.^ N Zn (21) 

N 

Then a useful simplification occurs because the sums over the m^ac/.'s are now 
independent; moreover the factors er^ can be reconstructed by multiplying every 
number T diag , by e^™*^ , so that: 

Z 9-c = II eX P t eX P (P^diag.) X f diag . X T diag ] (22) 
diag. 

We therefore obtain the grand potential (multiplied by —f3) in the form: 

LogZ g . c . = e ^ nMa9 - X fdiag. >< ^diag. (23) 

diag. 

This is an exact formula, which gives the value of the pressure of the system 
(multiplied by its volume and divided by the temperature). 

The second summation which can now be done consists in grouping together 
the contributions of all U-C diagrams that have the same "frame" (or "skeleton); 
in other words we sum the diagrams which have the same topology and differ 
only, inside the horizontal lines that represent the permutation sycles, by the 
lengths of the intermediate chains of U\s that connect together the particles 
contained in the C/j's with I > 1. Beside the fact that this operation turns out 
to be mathematically simple - it merely leads to the introduction of fractions 
of the operator U\ as we will see in the next section (see also the discussion 



of §2 of ref. [10]) -, it is also indispensable from a physical point of view: we 
have to make a summation over all lengths of intermediate horizontal chains of 
U\& in order to take into account an arbitrary degree of degeneracy. We call 
"classes" these groups of topologically equivalent U-C diagrams; classes may also 
be represented by diagrams (H-diagrams) which are, in a sense, simpler than 
the original diagrams since any indication of the length of the cycles has been 
removed. Examples are shown in figure 3 where dashed lines emphasize that a 
summation over cycle length is implied. We call 5 c ; ass the contribution of a class: 

^class ^ "j 6 ^ dla 9- x f d i a g X ^ diag. (24) 

diag. £ class 
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In terms of classes, ( p3[ ) becomes: 

Log Z g _ c , = -class (25) 

classes 

As (0), this is an exact formula, containing extensive quantities in both 
sides, and thus well adapted to approximations (as opposed to Z g , c ). We will see 
below that the first term of the summation gives the grand potential^of an ideal, 
(degenerate), gas, which we shall note "^ideai- The generalized Beth Uhlenbeck 
formula is contained in the second and the third term in the summation, which 
we shall note direct an d ^ moreover, in each of these two classes, one 

limits the summation of (|24]) to its first term (lowest order in Ui), one recovers 
the usual formula, valid only for weakly degenerate gases. 



5 Dilute degenerate systems 
5.1 Ideal gas 

We first check that the first class of diagrams reconstructs the grand potential 
(multiplied by —(3) of the ideal gas. This class, symbolized in figure 3-a, corre- 
sponds to the summation of the contribution of pure exchange cycles containing 
only C/i's, summed over any length k ranging from one to infinity. We know 
from (|18|) the counting factor, so that we just need to calculate the numerical 
contribution of every cycle. The result, proved below, is simple: 

r fc = ^ +1 Tr{[^f} (26) 

where: 

{1 for bosons 
— 1 for fermions 

This is again a consequence of the fact that the trace over the iV particles con- 
tained in the diagram is taken in a space that is simply the tensor product of k 
single particle spaces of state. Because the numbering of the particles does not 
affect the value of (it just changes the names of dummy variables), we can for 
convenience renumber the relevant particles from 1 to k. The effect of Ck is then 
to move particle 1 into the place initially occupied by particle 2, particle 2 into 
the place occupied by particle 3, and so on until one comes back to the place of 
particle 1. Introducing a complete set of states {| ip n >} in the one particle space 

9 More precisely, the logarithms of partition functions give the value of the grand potential 
multiplied by — (3. 
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of states^, one can then write: 
T fc = V k+1 E < 1 : Vnx I tfi(l) I 1 :( Pn 2 >< 2 : (f n2 I U x {2) \ 2 : <p n3 > x 



ni,n 2 ,...n k 



x....x < fc: Vnk | C7i(A:) | fc : <p ni > = V k+1 ^ {[Uif} 

(28) 

(the factor 7] k+l is equal to the parity of the cycle which enters in the definition 
of A for fermions). We now have to make the summation: 



:xj -y oo -y 



-"ideal £ 4 i~ - « £ ^ i 

k=l K k=l K 



which contains a well known series: 



E = E ^V +1 Tr {[U,]"} (29) 



9 S 4 

x + vy + y + + "' = ~ r?Log t 1 ~ ^ ^ 

We therefore get for the grand potential (multiplied by —0) of the ideal gas: 

Z id eai = -vTr {Log [l - ne?" U 1 ] } (31) 

This is the classical result. For instance, we can assume that Hq, the one-particle 
hamiltonian, is equal to P 2 /2m (kinetic energy of a particle in a box); by replacing 
in ( pip the trace by a sum over d 3 k, and U\ by its diagonal element e~^ n k ^ 2m , 
one immediately recovers usual formulas that are found in most textbooks on 
statistical mechanics. Indeed, the method that we have used is more indirect than 
the traditional method, but it gives a physical interpretation to the term in [Ui] k 
that is obtained by expanding the logarithmic function of ([H]): it corresponds to 
the contribution of all possible cyclic exchanges of k particles in the system. 



5.2 Second virial correction 

What happens now if we add the two following terms in (|25|), which contain one 
single U2 and no Ursell operator of higher order ? Let us start with the first class 
of diagrams, shown in fig. 3-b where the two particles in the U2 operator belong 
to two different exchange cycles, and which we will call direct diagrams. In the 
first diagram of this class, the two cycles are of length k = 1 (identities) and only 
two particles, unaffected by exchange, are involved; this simply introduces the 
contribution: 

rii = Tr 1)2 {C/ 2 (l,2)} (32) 

10 If the particles have internal states, the index n symbolizes at the same time the orbital 
quantum numbers as well as those characterizing the internal state. For instance, it the particles 
have spin /, a summation written as ^ n contains in fact two summations, one over orbital 
quantum numbers, and a second over (21 + 1) spin states. 
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The next diagram in this class corresponds to three clustered particles, two con- 
tained in the same U 2 and two in one permutation operator C 2 . The numerical 
value of this second diagram is: 



r 2 ,i = r 7 Tr lj2)3 {^2(l,2)?7i(3)C 2 (l,3)C 1 (2)} 
= r)Tn, 2)3 {U 2 (l, 2)C/ 1 (3)P ex .(l, 3)} 



(33) 



or: 



r 2 ,i = v E 

ni,n 2 ,n 3 

< 1 : W |< 2 : (f n2 | U 2 (l, 2) | 1 : ip n3 >| 2 : ip n2 >< <p n3 \ U x \ ip ni > 

(34) 

which provides the following reduced value: 

r a , 1 =7 ? Tr li2 {^(1,2)^(1)} (35) 

Similarly, one would calculate a contribution r 12 arising from the exchange of 
particles 2 and 3, and obtained by replacing in (|35|) U\(l) by Ui(2). More gener- 
ally, when U 2 clusters together k% particles, belonging to the same permutation 
cycle of length k\, with k 2 particles belonging to another cycle of length k 2 , the 
calculation of the effect of each of these cycles remains very similar to that of 
section 5Tj: now we have two particles that separately exchange with others, but 
the algebra of operators remains the same for each of them. We therefore get the 
reduced value: 



r^-y^Tn,. {U 2 (l, 2) pTi(l)]* 1 - 1 Wm^ 1 } (36) 



For this class of diagrams, according to (0) the counting factor / is simply 
1/2. The last step is to make a summation over all possible values of k\ and k 2 
after inserting an exponential of (3 times the chemical potential multiplied by the 
number of particles contained in the diagram: 



'—'direct ^ 

kito 



\ E e^ +k2) T klM (37) 



This operation can be done by using the relation: 

OO 1 

E M*- 1 = — (38) 
fc=l 1 '/ x 

One therefore introduces fractions of the U\S operators, which results in the 
expression: 



=1 



1 f e 13 ^ 



-direct ~ ^Tr 1>2 Wl,2) ^ (39) 
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For the second class of diagrams, exchange diagrams shown in figure 3-c, the 
two particles contained in U 2 are intermixed inside the same circular permutation. 
The first exchange diagram corresponds to the two particles contained in the same 
transposition: 

r£ = Tr 1>2 {U 2 (l, 2)C 2 (1, 2)} = Tr li2 {U 2 (l, 2)P CX .} (40) 
The second will contain three particles: 

r^ 2 = Tr w {U 2 (l, 2)L7 1 (3)<7 3 (1, 2, 3)} (41) 

which is equal to: 

F l*2 = E < 1 : l< 2 : <^ r 



m,n 2 ,n.3 

[/ 2 (1, 2) | 1 : ip n2 >| 2 : y2 n3 >< ip nz \U 1 \l: <p ni > 



(42) 



Now, we can use the equality: 

< 1 : ip ni |< 2 : ip n2 | U 2 (l, 2) \ 1 : ip n2 >\ 2 : ip n3 >= , . 

< 1 : <^ m |< 2 : ip n2 | C/ 2 (l, 2)P CX . | 1 : (p m >| 2 : (p n , 2 > { ' 

which allows us to get the same summation over indices as in equation ( |34"D and 
to obtain the reduced value: 

r^ 2 = Tr 1)2 {[/ 2 (l,2)P ex A(l)} (44) 

Another, very similar, term occurs if the circular permutation 6*3(1, 2, 3) of 
( |4"TD is replaced by C 3 (l, 3, 2); the calculation can easily be repeated and provides 
the result: 

r"i =Tr 1)2 {U 2 {l,2)P ex U 1 {2)} (45) 

From the preceding equations it is not difficult to see that the generic term 
of this second class of diagrams is obtained from ( p6|) by a simple replacement 
of U 2 by the product U 2 P ex .. Inserting a P ex . into (|39"D therefore provides 5g Xc/l 
Finally, the value of the grand potential (multiplied by —/3), to first order in U 2 , 
is given by: 

Lo g z, c . = + r ri , 2 {^(1, z) ll + f*\_J PUi(1)1 _, le l Ui{2) } 

(46) 

This result is valid within an approximation which is basically a second virial 
treatment of the interactions, while it contains all statistical corrections. The 
formula remains therefore valid if the degree of degeneracy of the gas is significant. 
Nevertheless, as pointed out for instance in section 2.1 of ref. [E| and in ref. [13 



virial series (even summed to infinity) are no longer appropriate beyond values 
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where the density exceeds that of a phase transition; this is because singularities 
in the thermodynamic quantities occur at a transition (in the limit of infinite 
systems). Therefore, for bosons, the validity of ([16|) is limited to non condensed 
systems. The discussion of what happens when a Bose Einstein condensation 



takes place will be given in [|TT 



Equation ([4g) is a generalization of the Beth Uhlenbeck formula, and reduces 
to it if the two denominators containing UiS are replaced by one, an operation 
which is valid in the limit of low densities where is small. To see the equiv- 
alence between the low density limit of relation (46 ) with the usual value of the 



second virial correction, we start from the definition of the Wk functions given in 
equation (14.35) of reference ||, which for k = 2 becomes^: 

W 2 (r u r 3 ) = 2(A T ) 6 £ I *rjrr,r 2 ) | 2 (47) 

where the functions \I/ ns (ri,r 2 ) are a complete set of stationary states for the 
system of identical particles (they are properly symmetrized); the thermal wave- 
length is defined by: 

T y/2nmkBT ^ 
Relation (Wfh can be transformed into: 



W 2 ( ri ,r 2 ) = 2(A T ) 6 < 1 : n,2 : r 2 | £„ a | ^ ns X | e^ H | 1 : r ls 2 : r 2 > 

(49) 

Inside this equation appears a closure summation over the symmetrized states of 
the system so that: 

W 2 (ri,r 2 ) = 2(A T ) 6 <l:r 1 ,2:r 2 |^e-^|l:r 1 ,2:r 2 > (50) 

(for fermions, S is replaced by A). From this function ref. || defines the Ursell 
function U^(ri,r 2 ) by: 

U? (n,r 2 ) = W 2 (r 1 ,r 2 ) - ^( ri )^i(r 2 ) = AW 2 (t u t 2 ) + U^ 9 ™^^) (51) 

where AW / 2 (r 1 ,r 2 ) is the difference between the values of W^Cn.,^) w hh and 
without interaction potential, and U 2 dealgas {^i, r 2 ) the value of the Ursell function 
for a system of two free particles. Equation (14.49)0 of |§ shows that the second 
virial coefficient in the expansion of the grand potential is half of the integral of 
17^(^1, r 2 ) over the variables r x and r 2 , multiplied by (Ay)~ 3 and the inverse of 
the volume. But the contribution of U2 dealgas {ri, r 2 ) is automatically contained in 
LogZg c ., so that we can concentrate on AW / 2 (r 1 , r 2 ) only. Because the summation 

11 We assume for simplicity that the particles have no internal state (or, if they do, that they 
are all in the same internal state). 

12 In the second edition, this equation is numbered (10.49). 
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over ri and r 2 can be written as a trace in a space which is the tensor product 
of two one particle state spaces, and because: 



S,A= l -[l± P ex ] (52) 



we get the result: 

[1 + TjPex] 



J d 3 n d 3 r 2 AW 2 v (n,r 2 ) =2(A T ) 6 Tri, 2 j[/ 2 (l,2) 



(53) 



where £7 2 (1, 2) is the difference between the exponentials of the interacting particle 
hamiltonian minus that of free particles, which is precisely our definition (|3|). We 
therefore recover our result ([!§), provided the denominators l — rje^^Ui inside the 
trace are replaced by one. 



5.3 Discussion; reduced density operators 

The only difference between equation fl46|) and the Beth Uhlenbeck formula arises 
from the presence of the two fractions inside the trace. Since: 

1 _ e^U x 
1 - rje^Ux ~ 1 + ^l + rje^Ux ^ 

they are actually nothing but operatorial forms of the usual Fermi or Bose factors 
(1 + 7]f) that appear, for instance, in the collision term of the Uehling Uhlenbeck 
or Landau kinetic equation (with the usual notation / for the distribution func- 
tion) . In the present case, nevertheless, because i7 2 and U\ do not commute in 
general, the operatorial character of the fractions is relevant. For free (but mutu- 
ally interacting) particles in a box, the eigenvectors of U\ are plane wavesPl while 
those of U2 are different since they involve correlations between the particles. 

If the trace in the right hand side of (f5^) is calculated in the basis of plane 
waves, the correction is expressed as an integral containing the diagonal elements 
of U2 between such plane waves. The simplest situation occurs when all these 
diagonal elements have the same sign. Then, for bosons, because the eigenvalues 
of the fractions are larger than 1, the effect of degeneracy is always to enhance the 
effects of interactions; this is physically satisfying since the Bose Einstein statistics 
tends to favor situations where particles are close. In particular, if the potential 
is attractive and if there are two body bound states (molecules), their weight 
will be increased with respect to what is would be in the usual Beth Uhlenbeck 
formula. If, on the other hand, the diagonal elements of U2 between plane waves 
change sign when the relative momentum of the two particles changes, which may 

13 If the particles are subject to the effect of an external potential (atoms in a trap for 
instance), the eigenvectors of U\ are not simple plane waves, but the essence of our analysis 
remains valid. 
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happen if the potential has attractive as well as repulsive parts, more complicated 
cancellation effects may take place in both the usual Beth Uhlenbeck formula 
and its generalization, so that no general prediction on the sign of the effect of 
degeneracy is possible (except for the contribution of bound states which remains 
enhanced as above). 

For fermions, the eigenvalues of the fractions are between and 1 so that the 
effect of statistics are just the opposite of what they are for bosons: they tend to 
reduce the effects of the interactions, except if mutual cancellation effects take 
place. Moreover, if the system is strongly degenerate, the effect of the product of 
the two fractions is to cancel the contribution of all matrix elements corresponding 
to particles inside the Fermi sphere, leaving only interactions between particles 
near the surface (or outside) of the Fermi sphere. This applies to bound molecules, 
which introduce a contribution containing the scalar product of the bound state 
wave function by all plane waves outside of the Fermi sphere, exactly as in the 
Cooper problem. 

The applicability of pBj ) is not limited to the equation of state of the quantum 
system, it can also be applied to the evaluation of partial density operators, for 
instance the one particle density operator. If we setj^]: 

Z7,(l,2,..0 = [EMI) x Ut(2) x ... x U X {1)] 17(1, 2,. i) (55) 

we can choose U\ and the E/j's, instead of the E/j's, to characterize the correlations 
between the particles; by definition U\ is equal to one. Assume now that we keep 
all E/z's constant and vary U\ according to: 

d\J\ = de \ <f >< f | Ux (56) 

What is then the variation of Z? Inserting this relation into fl55]) , then into @ 
and finally into (|IJ), we get: 

f S i=N 1 
dZ N = de Tr h2 ,3,...N \ A K N J2 \ i : : V 9 I J ( 57 ) 

This expression is nothing but the product of Zn by the average value of the one 
particle operator Y%Zi I i '■ '■ V 9 I- Therefore, in terms of the one particle 
density operator^ pj, we obtain: 

Log Z N = NTr{\ ip) (if | pj} = N{ip \ Pl \ if) (58) 

de 

(the derivative with respect to e is meant at the value e = 0). We can therefore 
obtain reduced density matrices from this kind of operatorial variations. If the 

14 If one wishes to introduce hermitian operators only, one can symmetrize (^5|) by putting 
square roots of C/i's on both sides of U. 

15 Here we take the convention where the operator pj is normed to the number N of particles 
(instead of one). 
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grand canonical ensemble is used, the many-particle density operator is defined 
by: 

P 9 c = [Z a r 1 T.^ N { S A K N S A ) (59) 



N 

ac. 



and we find the following relation involving the one particle density operator pf : 

^ Log Z gc = Tr {pf | <p)(<p \} = (<p\ pT I <p) (60) 

This formula shows that pf c is the "operatorial derivative" of Z gc (we have only 
considered diagonal variations of U\ but, since the basis of states | ip) is arbitrary, 
a generalization to off-diagonal variations is easy). 



For the ideal gas, all Ui are zero if I > 2. We can then simply use (|31[) and take 
a derivative with respect to e by choosing in succession for | <p) all eigenvectors 
of pf ■ In this way we obtain: 

pr = i - n^u, (61) 

which is a well known result. When interactions are added, one should replace 
( |3lD by (H) with the substitution (|55D for [7 2 - One then gets a modified expression 
of pi where contributions arise from the derivatives of the denominators in ( [46]) as 
well as that of the numerator (through the variation of U2), that we do not write 
here for brevity. This shows how the method is adaptable to the evaluation of 
microscopic quantities; if one is interested in microscopic correlations, two body 
density operators can be obtained by a similar method. This is a specificity of 
our technique, as compared for instance to that based on the pseudopotential; 
the latter is built in order to reconstruct the eigenvalues of the hamiltonian, and 
hence the partition function, but not the wave functions, in particular certainly 
not at short relative distances. 

A final remark is related to the convergence of the power series that we have 
summed into fractions of the U\ operators (the remark applies for fermions only). 
As noted by Kahn and Uhlenbeck |J, when the chemical potential of a system of 
fermions becomes positive, the virial series diverge, and the equation of state is 
obtained by a continuation of a an analytic function. Here we observe the same 

phenomenon: the series in e^^Ui that we have summed over the size k of the 
exchange cycles becomes divergent when /1 > 0; nevertheless the sum remains a 
regular function and, for fermions, ([46]) has no singularity. 



6 Conclusion 

The basic idea of the method of Beth and Uhlenbeck, or more generally of Mayer 
cluster expansion methods |16|, is to reason in terms of functions which directly 
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give local approximations of the thermal equilibrium. Here we start from Ursell 
operators which give rise to various contributions in the form of U-C diagrams. 
This leads to equation (|25|) which provides an exact expression of the grand 
potential of the system, and expresses it as sum of various terms arising from 
H-diagrams, already containing a sum over all possible sizes of intermediate ex- 
change cycles. Each term is obtained as an integral (a trace) over a finite number 
of variables. The general expression is valid for dilute or dense systems as well, 
such as liquids or even solids, but for gases it becomes simpler because it can 
be truncated more abruptly. Indeed, the generalization of the Beth Uhlenbeck 
formula is obtained by limiting the summation to the first three terms only (the 
first corresponding to the ideal gas). 

In its spirit, the method that we have used is close to the "binary collision 



approximation of Lee and Yang [14| [15]" ; they also introduce the Ursell functions 
of an auxiliary system obeying Boltzmann statistics, and they establish relations 
between the Ursell functions of the two systems (see also ref. |17j for a discussion 



of this type of method). But these relations acquire a rapidly increasing com- 
plexity with the order of the Ursell functions, and since they do not use Ursell 
operators, it becomes necessary at some point to make approximations (integra- 
tions over the inverse temperature (3). In our method, we introduce operators 
which generalize Ursell functions, and we make no attempt to evaluate the Ursell 
functions of the system of identical particles. Rather, we calculate the complete 
contribution of each initial (unsymmetrized) Ursell operator to the symmetrized 
partition function; this explains why the final results are markedly different. 

Another point of comparison is the more recent calculation by Nozieres and 
Schmitt Rink JTSfl , who give a calculation of the thermodynamic potential for 
a system of fermions at low densities. In order to simplify the summation of 
the diagrams in their calculation, they assume that the matrix elements of the 
interaction potential are separable into a product of functions; moreover, for 
brevity, they explicitly include only one phase shift (s wave approximation) so 
that their result is less general than (HBI). Of course this does not mean that their 



method can not be generalized to fully recover (|46|), but we have not examined 
the question. In the same vein, a more detailed comparison between (f46|) with 
the results of Galitskii || for fermions and Belyaev |19[ for bosons, valid a zero 
temperature, would be useful. 
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APPENDIX I : GENERAL RULES FOR CONSTRUCTING DIAGRAMS 



Our convention for building the diagram associated to any particular term in 
the double summation (|9|) is the following: 

(i) we start from the Ui, or the Ui's, that are of highest order I — Im in this 
particular term, and identify the particle it (or they) contain(s) that has the 
lowest index number n min ; this particle is considered as the "base particle", and 
belongs to the "base cycle" from which all the rest of diagram will be drawn. By 
convention, the base particle is shown first (left position) and the base cycle is 
that at the lowest position in the figure; for instance, in figures 1-b and 1-d, this 
base is a three particle cycle. 

(ii) to continue in the construction of the diagram, we add a second generation 
of cycles. We first use the Ursell operator Ui M that contains the base particle and 
add the other permutation cycles which include the other particles in the same 
Ui M ; this is done in the order of increasing values for the index number of the 
particles contained in this Ui M , so that the order of the new cycles is clearly 
defined. This also defines, for each of them, a "secondary base particle" that is 
put first in the diagram^. 

(iii) we continue the addition of this second generation of branches in the 
diagram by moving along the base cycle and skipping all numbered particles that 
are in U\'s, until we reach one which belongs to an Ui with / > 2; we then add 
additional cycles containing the other particles inside this Ui. We use the same 
rule an in (ii) and we define "secondary bases" for the new cycles, so that their 
representation is also uniquely fixed. Going along all the base cycle in this way 
completes the first generation of additional cycles. 

(iv) Then we build in the same way the second generation, by starting in 
succession from all of cycles of the first generation in the order in which they 
were added, etc. until, eventually, the complete diagram is obtained. 

With these conventions, all cycles are individually identified, so that it makes 
sense to vary their lengths independently to generate all terms of (|9|); this will 
be useful below for the calculations which lead from the U-C to the H-diagrams. 
Figure 2 shows examples of diagrams; Appendix II gives more details on the 
conventions used in their representation (in particular, it is convenient to assume 
that, inside the trace, the series of Ursell operators are put before the series of 
cycles) . 

Now suppose that we reverse the question: starting from a given diagram, 
how do we identify its occurrence in the double sum? what is its weight in the 
summation that gives Log Z7 Assume that we throw randomly ndi ag . numbered 
particles into all available locations. With the precise rules that we have chosen 

16 If the connection is multiple, that is if the additional cycle contains several particles of the 
same base Ursell operator Ui m , the secondary base particle is that of lowest index number; all 
the other particles are then automatically located in the diagram by their order in the exchange 
cycle, and do not play a special role at this stage. 
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above, it is clear that double counting problems are avoided, but also that not 
all of these random configurations obtained are allowed. What is the proportion 
fdiag. of the configurations that are compatible with our conventions? Let us note 
Pi M the total number of operators Ui M (those of largest order Im) appearing in 
this particular diagram. The reasoning is as follows: 

(i) first there is a probability (pi M x Im)~~ that the right base particle will be 
obtained 

(ii) second, if Im > 3, in the construction of successive generations of cycles, 
either more than one cycle is added from the connections of the same t/j, or there 
is a multiple connection towards the same cycle (or both). Figure 1-e gives one 
example of such a multiple connection. In all these cases, additional factors (g) 1 
are introduced which account for the correct ordering of numbering of secondary 
base particles. 

Thus we obtain: 

fdiag. = [Pi M x l M x g}' 1 (62) 

(this formula is also valid if Im = 1, in which case pi M is nothing but the size 
k of the linear exchange cycle). For the generalization of the Beth Uhlenbeck 
formula, case (ii) never happens and the factors g's do not play any role; they 
would nevertheless enter the calculations in general, for instance in the discussion 



of Bose Einstein condensation of groups of three particles in [|TT 

APPENDIX II : TWO EXAMPLES 

The diagram shown in figure 2-a corresponds by definition to the following 
trace: 

7V 1)2 , 3j 4, 5 {U 2 (l, 2)U 2 (3, 5)C/ 1 (4)C 1 (1)C 3 (2, 3, 4)d(5)} (63) 

Our convention is that the Ursell operators are always put before the cycles; the 
notation C 3 (2, 3, 4) refers to a cycle where the particle numbered 2 replaces that 
numbered 3, that numbered 3 replaces that numbered 4, and that numbered 4 
that numbered 2 (the C\ 's do not produce any change in the positions of the 
particles). Equation (jB3|) gives what we call the "explicit value" of this particular 
diagram, but it can also be simplified into a "reduced value". This can be done 
by inserting into ( |63"D closure relationships and using summations to introduce 
products of operators whenever possible. In this case this leads to the expression: 

T n ^{U 2 (l, 2)^(2,5)^(2)} (64) 

In the explicit value, any numbered particle appears once and only once in every 
Ursell operator and every cycle; in the reduced form this is not necessarily the 
case. 

In a similar way the diagram shown in figure 2-b is defined by the explicit 
expression: 

Tr li2 , 3 ,4,5,6 {U 2 (l, 2)Z7i(3)Z7 2 (4, 5)f/ 1 (6)C 1 (l)C 3 (2, 3, 4)^(5, 6)} (65) 
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while its reduced value is: 

vTr lt2 , 5 {U 2 {1, 2)C/ 1 (2)C/ 2 (2, 5)^(5)} 



(66) 



FIGURE CAPTIONS 

Fig. 1: Examples of U-C diagrams. For an ideal gas, only linear U-C diagrams 
containing chains of one particle Ursell operators U\ occur, as shown in (a); a 
summation of the contribution of these diagrams over the length of the chain gives 
the grand potential (multiplied by —(3). The generalization of the Beth Uhlenbeck 
formula arises for the diagrams shown in (b) and (c) containing one single two 
body Ursell operator U2] for each of them, a summation over the lengths of the U\ 
chains is also necessary. Figure (d) shows an example of a diagram containing a 
three body U3 operator and three U% chains, with a counting factor g = 2. Figure 
(e) shows an example of multiple connection through an U3 operator, which also 
introduces a factor g = 2 in the weight of the diagram. For more details on the 
definition and counting of the diagrams, see Appendix I. 

Fig. 2: Other examples of more complex U-C diagrams. As in figure 1, hori- 
zontal lines correspond to exchange cycles containing U\ operators; in addition, 
vertical double lines symbolize U2 operators (triple lines would be used for U3, 
etc.). Formula fll9"D gives the weight of these diagrams; a more explicit calculation 
of these two terms is given in Appendix II. 

Fig. 3: H diagrams introduced by the summation of the diagrams of figure 1 
over the lengths of the U\ chains, according to formula ([24]). The dashed lines, 
which symbolize these summations, can be replaced by intermediate operators 
given by fractions 1/ [1 — tjzUi], where z = exp (3/j,; every operator C/j, with I > 2, 
remains explicit and, moreover, introduces a factor z l ; finally, the weights fdiag. 
must also be inserted in the value of the H diagram. 
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